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. Baryon Acoustic Oscillations (BAO) provide an important standard ruler which can be used to 

■ probe the recent expansion history of our universe. We show how a simple extension of the Om 

diagnostic, which we call Om3, can combine standard ruler information from BAO with standard 
candle information from type la supernovae (SNIa) to yield a powerful novel null diagnostic of the 
cosmological constant hypothesis. A unique feature of Om3 is that it requires minimal cosmological 
assumptions since its determination does not rely upon prior knowledge of either the current value of 
the matter density, flom and the Hubble constant Ho, or the distance to the last scattering surface. 
Observational uncertainties in these quantities therefore do not affect the reconstruction of Om3. 
We reconstruct Om3 using the Union 2.1 SNIa data set and BAO data from SDSS, WiggleZ and 
6dFGS. Our results are consistent with dark energy being the cosmological constant. We show how 
Om and Om3 can be used to obtain accurate model independent constraints on the properties of 
dark energy from future data sets such as BigBOSS. 
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on 

. ! 1. INTRODUCTION 

> 

The expansion of the universe appears to have undergone a dramatic change in its recent past with most observations 
t"**- ' suggestive of the fact that cosmic expansion began to accelerate when the scale factor a(t) of the Universe was about 
0.6 of its present value, i.e. at the redshift z ~ 0.7. Although cosmic acceleration is a dramatic phenomenon, evidence 
for it is indirect and stems from data observed along our past light cone and involving either standard candles (SNIa) 
lO ! or standard rulers (BAO, CMB) [H-i]. 

The raison d'etre behind cosmic acceleration remains unknown and the fact that it might signal the need for new 
physics is partly responsible for the high levels of activity marking this field. Possible drivers of acceleration range 
from the Einstein's early suggestion of the cosmological constant Q to more elaborate constructions collectively called 
Dark Energy (DE) which may equally well be based on the introduction of a new physical field (physical DE) or on 
modifying the Einstein General Relativity (GR) (geometrical DE), or both H,||. To distinguish between these widely 
varying alternatives, model independent diagnostic tools capable of differentiating different classes of DE models are 
of great help [1] . 

With this in mind, we introduced the Om diagnostic [l^, E| which successfully distinguishes evolving DE from the 
cosmological constant on the basis of observations of the expansion history H(z). The two-point Om diagnostic can 
be defined as follows 

Om(z 2 ; Zl )= ^"I'f^ h(z)=H(z)/H (1.1) 

(l + z 2 y - (l + zxf 

so that [13 

Om{z; 0) = Om(z) = fj z) '~ 1 . (1.2) 
(1 +z) 6 — 1 



For a spatially fiat universe, H (z) can be recovered from the luminosity distance via a single differentiation [12H1 
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H(z) 



dz \ 1 + z 



(1.3) 
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which makes H(z) more robustly determined than the equation of state 



w{x) 



(2x/3) d InH / dx - 1 

i - {H /Hyn am x* 



X = 1 + z 



(1.4) 



which involves knowing the second derivative oi Dl as well as the cosmological matter density fiom- Thus observational 
uncertainties in and flom muddy the reconstruction of w(z) to a much greater extent than H(z) and Om(z) [Tol[Tl| . 
The expansion history, H(z), can also be determined directly from ages of passively evolving galaxies [17|> the time 
drift of the cosmological redshift fill . radial BAO's (l9l - [2l| or from a combination of the Alcock-Paczynski test and 
galaxy redshift space distortions [22J. This opens up the possibility of determining Om(z) from a combination of data 
sets involving the low redshift universe (z < few), precisely where cosmic acceleration seems to have originated ! 
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FIG. 1: The Om3 diagnostic (|1.9[) is shown as a function of separation, zz — Z2, for z% = 0.2, 22 = 0.35 and Qo m = 0.27. This 
diagnostic is clearly able to discriminate quite well between evolving DE models and ACDM, since, for ACDM the value of 
Om3 stays pegged at unity, while for Quintessence/Phantom Om3 rapidly evolves to larger/smaller values as the separation 
zz — Z2 increases. 

The possibility of using Om as a null diagnostic follows from the fact that, for the cosmological constant 

Om(z 2 ; z\) = rW . (1.5) 

In other words, the value of Om is redshift independent for the cosmological constant, while for other models of cosmic 
acceleration Om(z) is redshift dependent . This important property of Om can be used to considerable advantage by 
defining the difference diagnostic 10] 

Om diS (z 1 , z 2 ) := Om(zi) - Om(z 2 ) (1.6) 

and the ratio diagnostic 

Om I3lt i {Zi, z 2 , z 3 , Zi) := — — -, . (1.7) 

Om(z i -.z 3 ) 

From (jl.5p one immediately finds that the following equalities must hold for the cosmological constant 

Omdi$(zi, z 2 ) = 0, Om ia tio{zi, Z2, z 3 , Z4) = 1 . (1.8) 
A departure of Om rat i from unity therefore serves as a 'smoking gun' or null test for concordance cosmology: ACDM. 
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If any two of the four redshifts in Om ra ti are identical then Om va ± lo reduces to the 3 point diagnostic OmZ 

fm rat i (zi, z 2 , zi,z 3 ) := OmZ{zi, z 2 , z 3 ) = — — r , (1.9) 

Om(z 3 ; z\) 



where Om'd = 1 for ACDM. 1 (See [27|, |29j, |30( for other discussions of the Om diagnostic.) 

In this paper we focus on Om3 and show that the main difference between Om3 and Om is that these two diagnostics 
test the cosmological constant hypothesis using different cosmological observables. Whereas Om is a powerful null 
diagnostic of the cosmological constant as it requires only a prior knowledge of the expansion history, h(z), one should 
note that deriving h(z) directly from cosmological observables in a purely model independent and non-parametric 
manner is not always a simple task. Although supernovae data do allow one to reconstruct h(z) in a model independent 



manner 



25 1 , deriving h(z) directly from large scale structure data seems to be considerably more difficult (see [27j for 
detailed discussions). In contrast to Om, the Om3 diagnostic is specifically tailored to be applied directly on baryon 
acoustic oscillation data, and depends on fewer cosmological observables than Om in this case. 

Combining BAO information from the SDSS, WiggleZ & 6dF surveys with SNIa information from the Union2.1 
data set, allows us to determine the three point diagnostic, Om3, at several independent redshift bins thereby placing 
robust constraints on the nature of dark energy. 



2. THE OM3 DIAGNOSTIC RECONSTRUCTED FROM SNIA AND BAO DATA 

A key role in the determination of Baryon Acoustic Oscillations is played by the 'dilation-scale' distance [33 



D v (z) 



D{zf 



cz 



H(z) 



1/3 



where D(z) 



Dl{z) 
1 + z 



cdz 
Wz) 



(2.1) 



and we have used Dl(z) = (1 + z) 2 Da(z) to relate the luminosity distance Dl to the angular size distance Da and 
also made the assumption that the universe is spatially flat. 

Other important parameters which can be extracted from BAO's include the acoustic parameter 



A(z) 



iooD v (z)vn Qm h 2 



cz 



(2.2) 



and the ratio 



d(z) 



D v {z) 



(2.3) 



where r s (zcMB) is the sound horizon at the epoch when CMB photons decouple from baryons. 
Equation (|2.1j) allows us to relate the expansion history, H(z), to Dy and D as follows 



H(z) 



czD(z) 2 
D v {zf 



(2.4) 



Independent measurements of Dy and Dl can therefore be used to reconstruct H(z), as discussed in [27] . 

In this paper we shall demonstrate that the ratio of the Hubble parameter at two redshifts can play a key role in 
cosmological reconstruction. To see this note that 



H{zi\ Zj) 



H{zj) 



~ Dy{zj) ~ 
_D v {zi)_ 



(2.5) 



[1] Note a passing similarity between the diagnostic pair {Om rat ; , Om^iff} and the Statefinder diagnostic [23(1 {r, s\, where r =a /aH 3 , 
s = (r — l)/3(q — 1/2). In both cases, for the ACDM model one finds {Om rat i , Om^is } = {bO} and {r, s} = {1,0}. However, unlike 
the statefinder pair, the values of Oralis & Om rat i are strongly correlated, as seen from 111.61 1 and 111.91 1. Additionally, neither of 
Om dif f , Om rat j is able to distinguish steady state cosmology (SS) (S1a = l,f2om = 0), from SCDM = 0,f!om = 1) while the 

Om diagnostic and the Statefinders are more successful on this score, since Om(z) = for SS while Om(z) = 1 for SCDM; similarly 
{r, s} = {1, 0} for SS while {r, s} = {1, 1} for SCDM (see also [H). 
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equivalently 

*<**>-(!) 

in other words, ratio's of BAO parameters Dy , A, d(z) are related to the ratio of the Hubble parameter. From (|2.6[) 
we find that H(zi\ Zj) does not depend either on Hq, \/flQ m h 2 , or even the CMB parameter J" s (zcmb). 

We therefore arrive at the following important result: the value of Dl together with the value of either of the BAO 
paramaters {Dy,A,d}, determined at three independent redshifts, is sufficient to evaluate the Om3 diagnostic and 
define a null test of the cosmological constant hypothesis ! 

This follows immediately from & (|1.9|) , since, from the definition of Omi in (|1.9j) it follows that 

Om3(z i; z 2 ; z 3 ) = 3 3 / 3 _ 3 , where ir = 1 + z, (2.7) 

where H(zi\Zj) is determined using (|2.5[) or (I2.6[) . We should note that interchanging redshifts allows us to define 
different variants of Om3 (zi, z 2 and z 3 should not be necessarily sorted in redshift). Large galaxy redshift survey's 
including SDSS, WiggleZ and 6dFGS have determined the BAO parameters at several distinct redshift points [J-Q- 
In this paper we shall utilize these results to reconstruct Om3, thereby obtaining useful constraints on the nature 
of dark energy. It should be stressed that the model independent method which we advocate does not rely on a 
knowledge of either the matter density ttom or the Hubble parameter Hq. Its scope and prowess is therefore likely to 
dramatically improve as new data is added to the burgeoning BAO inventory by dedicated upcoming surveys such as 
BOSS 0, BigBOSS [azl, J-PAS [33 and Euclid 







A( Zj ) 




_ z t 


D{zi) 
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A{z % ) 


Z 3 






_d(zj)_ 



(2.6) 



A. Cosmological Data Sets 



As we stressed above, in order to determine Om3 one only needs to know the luminosity distance, Dr 1 (zi), and 
any one of the BAO parameters Dy(zi) i A(zi), d(zi). In this paper we shall use the Union 2.1 supernova data set @ 
to reconstruct D L , and the SDSS DR7 Q, WiggleZ 0,0] and 6dFGS [f| determinations of the BAO parameters to 
determine OmZ. Below we briefly describe these data sets and our method of extracting Om3 from them. 

1. Union 2.1 SNIa data. The Union 2.1 set 0] consists of 580 type la supernovae sampling the redshift range 
0.015 < z < 1.414. This data set includes 175 SNIa at low redshifts z < 0.1. 

Union 2.1 supernovae allow us to determine the luminosity distance Dl(z) in a model independent manner 
following the efficient smoothing ansatz proposed in (25|. Namely, a smoothed value for the luminosity distance, 
D£(z), is constructed from the fluctuating 'raw' value implied by data, Dl(z), by smoothing the latter using a 
low pass filter F having a smoothing scale A 



Df 



D L (z')F(\z-z'\;A)dz' , 



(2.8) 



where F(\z — z'\; A) is a variant of the Gaussian filter Fq oc exp(- 
iterative procedure [25j 



2 /2A 2 ). Specifically, one follows the 



In Dl(z, A) = In D 9 L (z) + N(z) V [In D L { Zi ) - In £>£(>;)] X exp 



hr 



2A 2 



(2.9) 



N(zy 



E 



exp 



In" 



l + Zj 

l+z 



2A 2 



where D 9 L (z) is a 'guessed' background model which is subtracted from the data before smoothing, thereby 
ensuring that it is the noise that is smoothed and not the luminosity distance ! The background model is 
upgraded at each iteration step via D g L (z) — > D^(z) and it is found that convergence is reached within a few steps, 
and that cosmological reconstruction is quite insensitive to the initial guess value D 9 L (z), which is reassuring. 
One might note that cosmological reconstruction using the smoothing ansatz usually yields a better fit (improved 
X 2 ) when compared with reconstruction using parametric methods [27j . The sensitivity of this method can be 
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TABLE I: BAO distances, from [j| 



BAO sample 


z 


d(z) 


A(z) 


6dFGS 


0.106 


0.336 ± 0.015 


0.526 ±0.028 


SDSS 


0.2 


0.1905 dz 0.0061 


0.488 ±0.016 


SDSS 


0.35 


0.1097 ±0.0036 


0.484 ±0.016 


WiggleZ 


0.44 


0.0916 ± 0.0071 


0.474 ± 0.034 


WiggleZ 


0.6 


0.0726 ± 0.0034 


0.442 ±0.020 


WiggleZ 


0.73 


0.0592 ± 0.0032 


0.424 ±0.021 



further enhanced by making the reconstruction error-sensitive via the substitution [In D r, (zj ) — In D 9 L (zi)] — » 
[ln_Di,(zj) — In D 9 L (zi)} /cr]j L r z .) in (]2.9p ; see [25hf27j for more details, and 0, [27], HH for other applications of 
the smoothing method. 

2. Data from Baryon Acoustic Oscillations 

Baryon acoustic oscillations - a relic of the pre-recombination universe - have been measured at 6 redshifts: 
z = 0.106, 0.2, 0.35, 0.44, 0.6 and 0.73. At these redshifts, the 'distilled BAO parameters' A(z) and d(z) have 
been determined to good accuracy, with the BAO detection itself being at the level of 2 — 3o\ Table 1, which 
has been reproduced from [5j, summarizes the BAO dataset. 

Not all of the BAO data in Table 1 is statistically independent. As pointed out in measurements belonging 
to the following redshift pairs are correlated: z = (0.2, 0.35), z = (0.44, 0.6), z = (0.6,0.73), the correlation 
coefficient being 0.337, 0.369 and 0.438, respectively. By using a BAO parameter associated with one of the two 
redshifts belonging to a correlated pair, one has 3 relatively independent parameters in all, which is precisely 
the correct number to determine Om3 ! 

In this paper we use d(z) results from Table 1 associated with the following redshift triplets to determine Om3: 
(0.106,0.35,0.6), (0.2,0.44,0.73) and (0.106,0.35,0.73). Our main result, therefore, consists of finding three 
independently obtained values for the Om3 diagnostic constructed from (|2.7|) & (|2.6p and using the BAO data 
in Table 1 jointly with Union 2.1 SNIa data. 

B. Results 

In figure [2] we show results for two variants of Om3 using three uncorrelated measurements of the WiggleZ baryon 
acoustic oscillation survey. In the left column the results for the {z\\ z^'-, z%) variant of Om3 are presented, while the 
right columns shows results for the (zi\z-j,\Z\) variant. The spatially flat ACDM model (OmZ = 1 for all variants) 
is in good accord with the data. However, the quality of the data is not yet good enough to allow us to confront 
different cosmological models with data in a precise and, at the same time, model independent way. Propagation of 
errors results in large error bars on Om3 making it difficult to distinguish DE models from each other. It therefore 
appears that Om3> has the potential to be used as a model independent future probe of cosmological models when the 
quality of data improves significantly. In figure [3] we show the expected value of Om3 determined using 500 simulated 
realizations of the BigBOSS experiment [37J. One clearly sees that Om3 carries the potential to strongly discriminate 
between rival cosmological models with least a-priori assumptions being made about the early/late universe. 

While Om3 can be very useful for testing cosmological models with minimal a-priori assumptions, having reliable 
information from the early and late universe can help us to determine h(z) and, knowing the latter, one can easily 
reconstruct the Om diagnostic. In figure 2] one realization of BigBOSS data together with expected future measure- 
ments of Hq with 2% uncertainty is used to determine Om(z). As we see from this figure, the Om diagnostic provides 
a powerful means to discriminate between rival DE models. 
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3. CONCLUSIONS 



In this paper we introduce a new null diagnostic customized for reconstructing the properties of dark energy from 
BAO data. Om3, as a 3 point diagnostic of dark energy, is closely related to the Om diagnostic and follows the same 
general principles. OmZ is designed in a way that can be applied directly to BAO and supernovae data in order to 
falsify concordance cosmology. This is done in a completely non-parametric way. The importance of Om3 lies in the 
fact that it does not rely on a knowledge of f2o m , Hq, or the distance to the last scattering surface, which permits Om3 
to falsify the ACDM model independently of these quantities. OmZ shares a common property with its cousin the Om 
diagnostic, namely, the values of both Om(z) and Om3(z) stay pegged (at unity for Om3(z) and f^om for Om{z)) in an 
expanding concordance cosmology (ACDM). This property of Om3 serves as a null test of the cosmological constant 
hypothesis, since, if observations do indicate that Om3 evolves with redshift, then this would imply that w ^ — 1 
for the equation of state of dark energy. We have shown that current BAO and supernovae data are in agreement 
with standard ACDM, however the uncertainties on Om3 are still pretty high. This is mainly due to the qual ity of 
available BAO data (current BAO data alone cannot put tight constraints on the cosmological quantities |28]) but 
future experiments, such as BigBOSS, can lead to much tighter determinations of Om3. Independence of Om3 from 
any a-priori assumptions of the early universe as well as values of H and fiom are important salient features of this 
null diagnostic of concordance cosmology (ACDM). 

Finally we would like to highlight the main differences between the two null diagnostics of ACDM - Om and Om3. 
Om is quite clearly a powerful null test of cosmological constant, and as it requires only a knowledge of h(z), it seems 
straightforward to use Om together with parametric or non-parameteric methods of reconstruction. However deriving 
h(z) directly from cosmological observables in a purely model independent and non-parameteric manner is not always 
an easy task. Though using supernovae data one can still derive h{z) in a model independent manner (25j . deriving 
h(z) directly from large scale structure data seems to be more difficult [27]. In this paper we have shown that Om3, 
in contrast to Om, is specifically tailored to be applied to baryon acoustic oscillation data (directly through BAO 
observables) and that, in this case, Om3 depends on a fewer number of cosmological observables than Om. 
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FIG. 2: Two variants of Om3 determined using uncorrelated BAO measurements at three redshifts in conjunction with 
the luminosity distance reconstructed from Union 2.1 SNIa data. The redshift combinations Om3(0. 106, 0.35, 0.60) & 
Om3(0. 35, 0.73, 0.106) shown in the top left and bottom right panels have the smallest error bars and therefore provide the 
tightest constraints on DE from current BAO+SNIa data. The large errorbars in the middle panels for Om3 are mainly due 
to the poor estimation of d(z = 0.44) (see Table 1). The spatially flat ACDM model is in good agreement with the data. Note 
that these results do not rely on a knowledge of Ho m , Ho, or the distance to the last scattering surface. 
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FIG. 3: Two variants of Om3 derived using simulated realizations of the BigBOSS experiment assuming a fiducial ACDM 
cosmology. Horizontal lines represent different dark energy models with (top-down) w = —0.8, to = — f .0, w = —1.2. flom = 0.27 
is assumed for all models. Note that the determination of Om3 requires minimal cosmological assumptions since one does not 
require a background model to estimate Qom, Ho, or the distance to the last scattering surface. Note also the difference in 
scale of the y-axis relative to figure [2] 
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FIG. 4: The Om diagnostic is reconstructed for a single simulated realization of the planned BigBOSS experiment assuming 
a fiducial ACDM cosmology. A determination of h(z), and hence Om(z), from future BAO experiments can clearly help 
distinguish between rival models of dark energy. Note that this determination is based on values of Ho with 2% uncertainty 
expected by the time BigBOSS becomes operational. 



